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ABSTRACT 

, We investigate the hydrodynamic evolution of Gamma-Ray Burst (GRB) remnants 

at late times, 3> 1 yr, by extending the standard fireball model into the nonrelativistic 

Sedov-Taylor regime. We calculate the associated synchrotron luminosity as a 

function of remnant age. The radio luminosity of old remnants depends strongly 

qq ' on the power-law index p = dlogN/dloge of the kinetic energy distribution of 

shock-accelerated electrons. In addition, the detection probability of remnants in 

nearby galaxies depends on the GRB rate in the local universe. We predict that if 

p 2.6 and the BATSE-calibrated GRB rate per comoving volume is independent 

of redshift (despite the strong evolution in the cosmic star formation rate) then the 

VLA can discover a ~ 10 years-old GRB remnant in the Virgo cluster. Nearby GRB 

remnants could be searched for and resolved at radio frequencies, and then identified 

as distinct from supernova remnants based on the detection of the extended emission 

of high-ionization lines in follow-up optical observations. 
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^ ■ 1. Introduction 



Since the detection of redshifted spectral features in gamma-ray burst (GRB) afterglows, 
it has become evident that GRB events do indeed occur at cosmological distances (Metzger et 
al. 1997; Djorgovski et al. 1998; Kulkarni et al. 1998a; Bloom et al. 1998; Bloom et al. 1999). 
The enormous energy release E ~ 10 51 — 10 54 ergs implied by the observed fluences and the 
cosmological distance scale could lead to a relativistically-expanding fireball, which produces the 
prompt (<^ 100 sec) 7-ray emission due to collisions between its internal shells (Paczyhski & Xu 
1994; Rees k Meszaros 1994; Pilla & Loeb 1997; Kobayashi, Piran, & Sari 1997). The fireball then 
enters the afterglow phase at later times (^ 1 day) when the expanding wind decelerates due to its 
interaction with the surrounding medium (e.g. Meszaros & Rees 1997; Waxman 1997a, b). During 
the early afterglow phase, the blastwave is still ultra-relativistic, and is hence well described by the 
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self-similar Blandford-McKee solution (Blandford & McKee 1976). Within a month the expansion 
is no longer ultra-relativistic. After about a year the expansion becomes nonrelativistic and the 
gas dynamics is well approximated by the Sedov-Taylor self-similar solution (Taylor 1950; Sedov 
1959), similarly to a supernova remnant (SNR). Given the BATSE-calibrated rate of occurrence of 
GRBs in the universe (Wijers et al. 1998) and the lifetime of their remnants, one infers that there 
should be hydrodynamic fossils of GRB remnants in any spiral galaxy at any given time (Loeb & 
Perna 1998). In fact, a subset of the so-called HI supershells observed in the Milky Way and other 
nearby galaxies might be old GRB remnants (Efremov et al. 1998; Loeb & Perna 1998). Rhode et 
al. (1999) have identified some HI holes in the nearby galaxy Holmberg II that do not have optical 
counterparts as expected in alternative models involving multiple SNe or starbursts with a normal 
stellar mass function. 

There are two important physical differences which could in principle lead to ways of 
observationally distinguishing GRB remnants from SNRs. First, the radiation energy emitted 
by an isotropic GRB explosion is estimated to be up to 4 orders of magnitude higher than that 
released in a supernova (Kumar 1999). Hence if a nearby expanding remnant is discovered and 
the explosion energy (as derived from the expansion speed and the Sedov-Taylor solution) is too 
great to be attributable to a single supernova, then one must conclude that we are seeing the 
remnant of either a GRB or multiple supernovae. But a GRB releases its total energy promptly, 
while the explosion time of multiple supernovae in a star-forming region is not expected to be 
synchronized to better than ~ 10 6 years. Hence, GRB remnants should be much more energetic 
than SN remnants at early times. The existence of energetic X-ray remnants was recently 
inferred in deep ROSAT X-ray images of M101 (Wang 1999); but future observations with the 
XMM and Chandra X-ray satellites are necessary in order to identify the nature of these sources 
conclusively through high-resolution imaging and spectroscopy. Second, GRB afterglows include 
a strong UV flash, which on the timescale of a few hundreds of years creates an ionized bubble 
of radius ~ 100 n\ pc, where n\ is the ambient density in units of 1 cm~ 3 (Perna & Loeb 
1998). At ages <J 10 4 years, the non-relativistic shock wave acquires a much smaller radius 
than the ionized bubble. This provides a potentially unique signature of young GRB remnants; 
namely, a relatively compact expanding blastwave embedded in a much larger ionized sphere. The 
optical-UV recombination lines from highly-ionized species in this region should provide a clear 
discriminant for recognizing GRB remnants (Perna, Raymond h Loeb 1999). The embedded blast 
wave emits both thermal and non-thermal radiation; the latter being an extension of the afterglow 
due to synchrotron emission by shock-accelerated electrons. 

In this Letter we calculate the synchrotron radio emission from the expanding GRB shock 
in old GRB remnants. Our goal is to examine whether a survey in the radio could efficiently 
identify old GRB remnants. Once a candidate has been identified in the radio, follow-up optical 
observations could search for the special recombination lines which are generic of GRB remnants. 
In §2 we review the blastwave hydrodynamics as extrapolated into the Sedov-Taylor regime, and 
calculate the synchrotron flux as a function of age, as well as the expected source counts in the 



-3- 



Virgo cluster. Finally, we discuss our conclusions in §3. 



2. Hydrodynamics and Emission 

While the hydrodynamic evolution of GRB remnants through the early afterglow phase is 
well-described by the ultra-relativistic Blandford-McKee (1976) solution, at times much greater 
than a year we may to a good approximation use the non-relativistic Sedov- Taylor solution. It 
has been shown (Huang, Dai, & Lu 1998) that the two regimes may be smoothly matched to 
one another. We assume that the kinetic energy remaining in the fireball is comparabale to the 
total energy Eq which is radiated away in 7-rays, and that the energy release is isotropic; the 
effects of beaming in the initial energy release will be discussed in §3. The mean atomic weight of 
the surrounding (interstellar) medium is taken to be lAm p , where m p is the proton mass. The 
appropriate scaling laws for the shock radius r s and velocity v s in the non-relativistic regime are 
then 

r s = 1.17 (y!^) ' 2 = 1-56 x 10" nr°- 2 ^ # cm, (1) 

v s = 0.47 {^-^j' 2 = 1.99 x 10 10 n r °- 2 t£* cm s" 1 , (2) 

where n = 1 ri\ cm~ 3 is the pre-shock particle number density, Eq = 10 52 -E52 erg is the total 
burst energy, and t yr is the age of the remnant in years. 

The state of the gas behind the shock is defined by the Rankine-Hugoniot jump conditions. 
For a strong shock, the post-shock (primed) particle density, bulk velocity, and energy density are, 
respectively, 

ri = An, (3) 
v' = \v s , (4) 

u' = ^nm p v 2 s , (5) 

where all quantities are measured in the frame of the unshocked material. 

As usual, we assume that the electrons are injected with a power-law distribution of kinetic 
energies, 

/(e) = (p - 1) e-f , e > e m , (6) 

where f(e)de is the fraction of shock-accelerated electrons with kinetic energies in the range 
(e, e + de). We assume a single power-law index at all energies, for simplicity. The spectral break 
that might exist at the transition to nonrelativistic electron energies occurs well below the electron 
energy which is responsible for the relevant radio emission at ~ 1 GHz. 
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The post-shock magnetic and electron energy densities are assumed to be given respectively 
as fixed fractions and £ e of the total post-shock energy density: 

B' 2 

v! = n 'e m = £ e u , (8) 

p-2 

where m e is the electron mass. The magnetic field strength and minimum electron energy are thus 
given by 

r>i _ n i a c0.5 „0.3 z, 

P ~ 2 c ^-0.4 ^0.4 + -1.2 
p—1 



B' = 0.14 e b 5 < 3 EH t;r G, (9) 
117 ^ £ e nr a4 £ 5 ° 2 4 i- 1 " 2 MeV. (10) 



Before calculating the flux, we need to determine whether the Sedov-Taylor remnant is 
optically thin to synchrotron self-absorption. The shocked region has a thickness r]r s , with 77 ~ 0.1 
from particle number conservation. We will adopt the value r] = 1/15, which is consistent with the 
density profile of the Sedov (1959) solution, and provides a number of shocked electrons which is 
a fair fraction of the total number swept up by the blast wave. The self-absorption optical depth 
at a photon frequency v is given by (cf. Rybicki & Lightman 1979, p. 190) 



y/3e*nB ± rir 8 fp 11\ fp 1 W e m V / JSeBj_ \^ 

Here, m e is the electron mass, T(x) is the gamma function, and B± ~ B' is the component of the 
magnetic field perpendicular to the electron velocity. Substituting the numerical factors into (11) 
gives 

= ^1 f^)°- 5P a{p) e 0.5 + 0. 2 5 P n 1.5-0. 2 5 P ^ £-1.6, (12) 

^ 9 V ^9 / 

where v = 10 9 fg Hz, and 

e-^fr r (5 + TH5 + 5) 

is of order unity. Thus, for typical power-law index values 2.1 < p < 3.2 (Li & Chevalier 1999), 
we find that t{v) <C 1 for 1/ = 1 GHz at an age of ^> 1 yr. 

The flux from old remnants is therefore given by 

K _ x ^ (14) 

where -D is the distance to the remnant, and P v is the volume emissivity (in erg s _1 cnr 3 Hz -1 ), 
given by (Rybicki & Lightman 1979, p. 180) 

_ 4V3e 3 rLB ± p- 1 /p 19\ fp 1 \ / e m y 1 / 3e£? ± \Q-5(p-i) 

m e c 2 p+1 U 12 / U 12/ Vm e c 2 / \27rm e ci// ' 1 j 
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Note that we are justified in using the synchrotron formula since the electrons which emit at GHz 
frequencies are ultra-relativistic: 

(^sf )° 5 = 40 45 ^°' 25 nr °' 15 E ^ ^ (16) 



TJlpC 2 



Coulomb collisions thermalize the electrons only at non-relativistic energies, well below the 
energies of interest for this discussion. Hence, we are justified in using the power-law spectral 
shape in equation (p~J). The flux should eventually show a cooling break at high frequencies, due 
to the fact that electrons with more than a threshold energy e c will radiate away their energy 
faster than the dynamical time. The synchrotron cooling time is given by t c = 67r(m e c 2 ) 2 /<7Tcei? 2 _, 
where ot is the Thomson cross section. The cooling frequency is thus given by 

- (Stf - M * 1011 09 ^ - <»> 

For remnant ages <^ 10 7 yr, the cooling cutoff occurs well above 1 GHz. Thus, the emission at 
GHz frequencies should fall on the simple P v oc z/ 1_p )/ 2 part of the spectrum. 



We may now calculate the numerical value of the expected flux 

0.5(p-l) 

'Mpc 

where Dy\^ c is the distance to the source in Mpc, and 



F v = f^f 5 ^ m £*K14*) g-X ,0.95-0.25, ^3+0.5, ^1-1.5, Jy> (lg) 



*)-e^)e-EirK5 + s) r (5-5) 



is of order unity. 



In figure 1 we plot the flux at 1.6 GHz from equation (18) as a function of remnant age t yI 
for p = 2.1 (F v oc t -105 ) and p = 3.2 (F u oc t -2 ' 7 ), the two values which bracket the range of 
power-law indices seen in radio supernovae (Li & Chevalier 1999). We have assumed that the 
sources reside in the Virgo cluster, at a distance of -Dmpc = 16. We also assumed sub-equipartition 
energy density of the magnetic fields and nonthermal electrons in the post-shock gas with = 0.1, 
£ e = 0.1; a typical interstellar medium density, n\ = 1; and rj = 1/15. In panel (a), we take 
E52 = 0.7, derived from the BATSE catalog for a source population with a redshift-independent 
burst rate per comoving volume, while in panel (b) we use .E52 = 14, which applies for population 
which traces the star formation rate in galaxies. For comparison, the horizontal line shows the 
5cr sensitivity threshold, .Fvla = 70^Jy, for the Very Large Array (VLA) in the 1.6 GHz band, 
using a one-hour integration in the VLA's most extended configuration (L. Greenhill 1999, private 
communication). Clearly the maximum age for a detectable remnant, and hence the number of 
detectable remnants, depends strongly on p. The vertical line, at t = 600 yr and 10 5 yr in panels 
(a) and (b) respectively, indicates the age at which we expect to detect at least one remnant in 
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the Virgo cluster. The intermediate plots of F u (t) in panels (a) and (b) are for p = 2.6 and 2.3, 
the respective maximum electron slopes for which a 600 or 10 5 -year-old Virgo remnant would be 
detectable by the VLA. 

To calculate the expected source counts, we note that the local GRB rate per L± (where L± 
is the characteristic stellar luminosity per galaxy in the local universe), assuming no beaming, 
is estimated to be T » 2.5 x 10~ 8 L+ 1 yr -1 (Wijers et al 1998). This value is derived assuming 
that the GRB rate traces the star formation history of galaxies; for a non-evolving burst rate, the 
inferred value is 150 times higher. The number of remnants per L* younger than t yT [and hence 
brighter than F u (t yT )] at a given time is then Tt yr . The best place to search for optical emission 
from old GRB remnants is in the Virgo cluster (Perna, Raymond, & Loeb 1999). There are ~ 2500 
galaxies brighter than B = 19 in this cluster; at a distance of 16 Mpc, this limit corresponds to 
an absolute magnitude magnitude Mb = —12. For typical Schechter (1976) function parameters, 
a = —1 and M* = —19 (Loveday et al. 1992; Marzke et al. 1994), this yields a total luminosity 
for the Virgo cluster of Lyir = 430L*. Thus, for a non-evolving GRB population, we need to look 
for remnants as old as t = (T ■ Lyir)^ 1 = 600 yr to be reasonably confident of seeing at least 
one remnant in the Virgo cluster. For a GRB population which traces the cosmic star formation 
history, we need t = 10 5 yr. As illustrated in Figure 1, for the burst parameters we have assumed, 
remnants this old could be detected by the VLA only for an electron power-law index p < 2.6 in 
the non-evolving case and p < 2.3 in the evolving case. 

For deeper volume-limited searches, e.g. of volumes probed by the Sloan Digital Sky Survey 
(SDSS; Gunn & Weinberg 1995), the increase in the number of galaxies surveyed N oc D s 
dominates over the decrease in detectable ages t oc D 2 /( 2 - 1 ~ 1 - 5 p) ) obtained by solving equation (18) 
for t yr . Thus, although at a distance of 250 Mpc we need to look for remnants which are younger 
than 10 2 years, the galaxy luminosity within this volume is ~ 10 5 L*, leading to the prediction of 
at least one detectable young remnant for p = 3.2 and many more for lower p values. 



3. Discussion and Conclusions 

We have calculated the expected synchrotron flux from ^ 10 5 -year-old GRB remnants in the 
Virgo cluster. Although the most revealing signature of a GRB remnant may be the presence 
of optical-UV recombination lines from high-ionization states of metals (Perna, Raymond, & 
Loeb 1999), it might be easier to search for nearby GRB remnants in the radio due to the lower 
background noise (from the sky plus the host galaxy) at radio frequencies. For a non-evolving 
GRB population, we find that if the spectral index of the shock-accelerated electrons p < 2.6, then 
one could find at least one ~ 600-year-old remnant in the Virgo cluster at the VLA detection 
threshold. Such a remnant should be characterized by a synchrotron spectrum F v oc v~ ' 8 or 
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flatter. g 

Perhaps the most poorly constrained parameter of the GRB sources is the beaming fraction 
fb, the fraction of 4tt steradians into which the initial 7-ray emission is emitted; the actual event 
rate may then be enhanced to /fj~ 1 r. However, the total energy released in an event scales as foEo, 
and so the synchrotron flux from the remnant is proportional to y0- 3+0 - 5p . The decrease in the 
required energy due to beaming may be counteracted by the fact that the efficiency for producing 
7-rays in the initial event is very low (Kumar 1999). We note that even if the initial (impulsive) 
energy release is beamed, the deposited energy, foEo, will be isotropized at the onset of the 
non-relativistic expansion phase. In addition, for a uniform distribution of sources in Eucledean 
space, the number of detectable remnants younger than a given age should decline with decreasing 
f b as oc / 6 - 1 (/&°' 3+ °' 5p ) 3/2 = / 6 °- 75p ~ - 65 . 

Due to the as-yet uncertain nature of the progenitors, it is not clear whether we should 
expect the GRB rate to directly trace the star formation history. Hence we considered both the 
SFR-tracing and non-evolving cases in Figure 1. We summarize these results in Figure 2, which 
shows the upper bound on p as a function of explosion energy for the two cases. The factor of 150 
enhancement in the rate for a non-evolving population allows for a greater probability of seeing 
younger bursts with steeper spectral slopes. It is important to note that if a young (<; 10 3 yr) 
remnant is detected in the Virgo cluster, our results will strongly suggest that the local GRB rate 
is consistent with the non-evolving scenario, but not with the star-formation-tracing scenario. 
This is the only immediate way to constrain the local burst rate (short of monitoring the ~ 10 6 
SDSS galaxies over one year and searching for a GRB explosion). The young radio remnants 
detected in this case should be well-embedded in the bubble that was ionized by the initial UV 
flash, and should be detectable in optical recombination lines. If GRBs follow the global star 
formation history, however, we only expect to see a ~ 10 5 -year-old remnant in Virgo, which is 
large and not well-embedded, and will also not be easily detectable at optical wavelengths (Perna, 
Raymond, & Loeb 1999). 

Is it justified to assume that the only significant flux at 1.6 GHz comes from the freshly 
shocked electrons within a distance rjr s behind the shock? Might the electrons which were shocked 
at t ~ 1 yr contribute a substantial amount of flux when the remnant is 10 4 years old? The 
Sedov- Taylor similarity solution implies that the volume occupied by the material behind the 
shock front is increasing with time as r 3 . Hence, the relativistic energy densities of the electrons 
and magnetic fields scale adiabatically as u' b oc u' e oc r~ 4 . The flux at a fixed frequency v due to the 
old electrons is then F° ld oc t~°' 8p ° ld , where p Q id is the power-law index measured for afterglows 
at ~ 1 yr; typically p Q \d ~ 3. Using this value, a comparison with equation (|l^) yields a ratio 
poid^pnew _ t i.5 Pncw -4.5 j where ^new is the flux from new i y _ s hocked electrons. Thus, if p ncw > 3, 

the flux from the old electrons will dominate at times i > 1 yr. For values in the range observed 



3 Although it appears that the energy flux vF v oc v 2 diverges at high frequencies, recall (cf. Eq. [I?]) that there 
is a break in the spectrum at the cooling frequency u c , above which, F v tx v~ 3 . 
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in radio SNRs (Li & Chevalier 1999), we are justified in neglecting the contribution from the old 
electrons. 

Our extrapolation to the non-relativistic regime may be tested by continuing to monitor the 
event GRB 970508, which is at z = 0.835 and had a 1.4 GHz flux of 249 ± 60 /iJy at an age of 354 
days (Frail et al. 1999, in preparation). 

Our discussion assumed that GRB sources release their energy impulsively in the form of the 
ultra-relativistic wind that produces the early 7-ray and afterglow emission. It is possible that 
more energy is released in the form of non-relativistic ejecta that catches-up with the decelerating 
shock and re-energizes it at late times. In this case, if /& = 1 then the total hydrodynamic 
energy (and the associated synchrotron flux) may be much larger than we estimated based on the 
GRB energetics alone. However, the situation might be different if fj, <C 1. Recently, there have 
been several claims for a potential detection of supernova emission in the light curves of rapidly 
declining afterglows (Kulkarni et al. 1998b; Bloom et al. 1999; Reichart 1999). We emphasize 
that any association between supernova and GRB events could be explored more directly in the 
local Universe by examining the ionization and hydrodynamic structure of SN-like remnants in 
the interstellar medium of nearby galaxies. For example, one could search for extended ionization 
cones in young supernova remnants, as expected if the intense UV emission from the associated 
GRB afterglows is collimated. Complementary information about the shock structure and 
temperature can be obtained by observations in the radio band (probing the synchrotron emission) 
or the X-ray regime (probing thermal emission from the hot post-shock gas). In particular, 
follow-up observations of the energetic X-ray remnants discovered by Wang (1999) in M101 or the 
optically- faint HI holes discovered by Rhode et al. (1999) in Holmberg II, would be revealing as 
to the nature of these peculiar objects. 

The expanding spherical shock front will typically acquire a diameter of ~ 1" after 10 4 years at 
the distance of the Virgo cluster. With VLBI, one can achieve up to sub milli-arcsecond resolution 
and hence easily resolve the radio-emitting shock. Since the shocked gas occupies a thin shell, the 
source should be strongly limb-brightened and possibly highly polarized at the limb (Medvedev 
& Loeb 1999). In the simplest case of an isotropic point explosion, the radio-emitting region 
will appear embedded inside a much larger region which was ionized by the prompt UV emission 
from the GRB and which emits optical-UV recombination lines (Perna, Raymond, & Loeb 1999). 
This distinct structure of a shock embedded in an HII region with high ionization states of heavy 
elements, is unique to GRB remnants since only the optically-thin wind of a relativistic GRB 
fireball can give rise to the intense hard radiation which produces a highly-ionized bubble out to 
large distances ~ 100 pc, in front of the shock. This is to be contrasted with ordinary supernovae, 
in which the UV emission from the optically-thick envelope is suppressed above the thermal cutoff. 

We thank Lincoln Greenhill for useful discussions. This work was supported in part by NASA 
grants NAG 5-7039 and NAG 5-7768. 
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Fig. 1. — Observed remnant 1.6 GHz flux as a function of age (solid lines) at the Virgo cluster 
distance D = 16 Mpc, for (a) a non-evolving GRB population and (b) a star-formation-tracing 
GRB population. We show three different power-law indices; the upper and lower values correspond 
to the range of observed slopes in radio SNRs, while the middle value is such that one remnant 
should be visible in the Virgo cluster. Shown for comparison is the VLA 5a sensitivity level, 
-^Vla = 70 /uJy (horizontal dotted line), and the corresponding age limit [t = 600 yr in panel (a) 
and t = 10 5 yr in panel (b)], such that there should be at least one younger remnant in Virgo. We 
assume £5 = £ e = 0.1, and 77 = 1/15. 
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Fig. 2. — Maximum values of the electron power-law slope p, such that one remnant will be 
detectable in the Virgo cluster (assuming /& = 1). We show best-fit results for a non-evolving 
source population (upper left) and a population which traces the star formation history of galaxies 
(lower right). The widths of the horizontal line segments reflect the la uncertainty in the burst 
energy, from Wijers et al. (1998). 



